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ABSTRACT 

We present measurements of the small-scale (0.1 < r < 1 /i _1 Mpc) quasar two-point correlation function 
at z > 2.9, for a flux-limited (i < 21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The 
amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude 
is comparable to or lower than power-law extrapolations (with slope 7 = 2) from the large-scale correlation 
function of the i < 20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating 
quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Re- 
producing the large-scale clustering amplitude requires that the active fraction of the black holes in the central 
galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in 
satellite galaxies 0.1 < f s < 0.5 at z > 3. 

Subject headings: black hole physics - galaxies: active - cosmology: observations - large-scale structure of 
universe - quasars: general - surveys 



1. INTRODUCTION 

With the rapid progress in observational and computational 
cosmology in the past two decades due to dedicated sur- 
veys and numerical simulations, it is now possible to study 
the quasar population within the hierarchical structure forma- 
tion framework (e.g., Kauffmann & Haehnelt 2000; Volon- 
teri, Haardt & Madau 2003; Wyithe & Loeb 2003; Hopkins 
et al. 2008; Shankar et al. 2008, 2009; Shen 2009). If lumi- 
nous quasars are the progenitors of the most massive galaxies 
today, then they occupy the rare peaks in the initial density 
fluctuation field, i.e., they are biased tracers of the underly- 
ing matter distribution (e.g., Bardeen et al. 1986; Efstathiou 
& Rees 1988; Cole & Kaiser 1989; Djorgovski 1999; Djor- 
govski et al. 1999). The quasar two-point correlation func- 
tion has now been measured for large survey samples to un- 
precedented precision (e.g., Porciani et al. 2004; Croom et 
al. 2005; Myers et al. 2006, 2007a; Shen et al. 2007, 2008, 
2009; da Angela et al. 2008; Ross et al. 2009). These stud- 
ies suggest that quasars live in massive dark matter halos of 
Mh-Ao ^ a few x \0 12 IT 1 Mq\ their bias relative to the under- 
lying matter increases rapidly with redshift. However, such 
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studies are unable to probe the smallest scales (r < 1 /r'Mpc), 
where matter evolves nonlinearly and the distributions of 
quasars within dark matter halos start to play a role in de- 
termining their clustering properties. This is because fiber- 
fed multi-object spectroscopic surveys usually cannot observe 
two targets closer than the fiber collision scale ~ 1'. 

Hennawi et al. (2006) compiled a sample of close quasar 
binaries at z < 3 by spectroscopic follow-up observations 
of candidates selected from the Sloan Digital Sky Survey 
(SDSS; York et al. 2000) imaging data. Using this binary 
sample, they measured the correlation function down to scales 
as small as R pmp ~ 15kpc, where R prop is the transverse sep- 
aration in proper units; they confirmed and extended pre- 
vious tentative claims (e.g., Djorgovski 1991) that quasars 
exhibit excess clustering on small scales (most notably at 
Rp IO p < 40kpc) compared with the naive power-law extrapo- 
lation of the large-scale correlation function. This small-scale 
excess clustering was confirmed by Myers et al. (2007b, 2008) 
in a more homogeneous sample, albeit at a lower level of "ex- 
cess". 

The large-scale quasar correlation function has now been 
measured at high redshift (z > 3, Shen et al. 2007), where 
quasars cluster much more strongly than their low redshift 
counterparts. It is natural to extend the work of Hennawi 
et al. (2006) to study the small-scale quasar clustering at 
z > 3. However, such investigations are challenging for two 
reasons: first, the number density of the quasar population 
drops rapidly after the peak of quasar activity at z ~ 2-3 (e.g., 
Richards et al. 2006); second, quasar pairs on tens of kpc to 
1 Mpc scales are rare occurrences - only < 0.1% of quasars 
have a close quasar companion with comparable luminosity. 
Hence a large search volume is needed to build up the statis- 
tics. Hennawi et al. (2009, hereafter Paper I) have, for the first 
time, compiled such a binary quasar sample at z > 3, which 
we use here to study the clustering of quasars on small-scales. 

The rareness of close quasar pairs is not in direct contra- 
diction with the major merger scenario of quasar triggering, 
because the probability that two quasars are triggered and 
identified simultaneously during the early stage of a major 
merger (i.e., with separations on halo scales rather than on 
galactic scales) is low in theoretical models (e.g., Volonteri 
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et et al. 2003; Hopkins et al. 2008). However, even a hand- 
ful of close quasar pairs will contribute significantly to the 
small-scale clustering amplitude because the mean number 
density of quasars is so low that the expected number of ran- 
dom companions on such small scales is tiny. Also note that 
although quasar pairs with comparable luminosities are rare, 
there might be more fainter companions (i.e., low luminosity 
AGN or fainter quasars) around luminous quasars (e.g., Djor- 
govski et al. 2007), as expected from the hierarchical merger 
scenario. 

In this paper we measure the small-scale quasar clustering 
at z > 3 using a set of 15 quasar pairs in the sample of Paper I. 
We adopt the same cosmology as in Paper I, with Q m = 0.26, 
f^A = 0.74 and h = 0.7. Comoving units will be used unless 
otherwise specified, and we use subscript plop for proper units. 

2. THE SAMPLE 

Our parent sample is the high-redshift binary quasar cata- 
log presented in Paper I. This sample includes 27 quasar pairs 
with relative velocity |Av| < 2000kms _1 at 2.9 < z < 4.5, 
down to a limiting magnitude i < 21 after correcting for 
Galactic extinction, selected over 8142 deg 2 of the SDSS 
imaging footprint prior to DR6 (Adelman-McCarthy et al. 
2008). The detailed target selection criteria, completeness 
analysis, and follow-up spectroscopy can be found in Paper I. 
To construct our clustering subsample, we first exclude eight 
pairs that failed to pass the selection criteria (for which the 
completeness cannot been quantified) described in §2 of Pa- 
per I, leaving 19 pairs. Second, the follow-up spectroscopic 
observations are the most complete out to an angular sepa- 
ration 9 « 60", because those targets were assigned higher 
priority for follow-up spectroscopy, and therefore we restrict 
ourselves to pairs with angular separation 9 < 60"; this re- 
striction excludes one pair at z < 3.5 and three pairs at z > 3.5. 
Our final clustering subsample thus includes 15 pairs with 
seven pairs at z < 3.5 and eight pairs at z > 3.5, with pro- 
jected comoving separations /? ~ 0.1 — 1 /i _1 Mpc and proper 
separations R plop ~ a few tens to a few hundreds of kpc (see 
fig. 8 of Paper I). 

The sparseness of the sample requires different techniques 
for measuring the clustering strength, from the traditional 
binned w p statistic (e.g., Davis & Peebles 1983). Here we 
adopt the Maximum-Likelihood (ML) approach used in Shen 
et al. (2009), as described below. We report our ML esti- 
mates and statistical uncertainties of the small-scale clustering 
in 32.11 the systematic uncertainties are discussed in 32.21 To 
reduce the impact of the selection incompleteness at z ~ 3.5 
due to stellar contaminants (see Paper I), and to explore red- 
shift evolution, we measure the small-scale clustering in two 
redshiftbins: 2.9 < z < 3.5 (low-z) and 3.5 < z < 4.5 (high-z). 

2.1. Clustering Measurements 

Here we recast the ML approach of Shen et al. (2009). 
We choose a power-law model for the underlying correla- 
tion function: £(r) = (r/ro.ML)~ 7ML - We then compute the 
expected number of quasar pairs within a comoving cylindri- 
cal volume with projected radius R to R + dR and half-height 
AH = 20 /r'Mpc. This half-height is chosen to reflect our ve- 
locity constraint in defining a quasar pair and to minimize the 
effects of redshift distortions and errors. Assuming Poisson 
statistics, the likelihood function can be written as: 



(!) 



where // = 2irRh(R) dR is the expected number of pairs in the 
interval dR, the index i runs over all pairs in the sample and 
the index j runs over all the elements dR in which there are 
no pairs. The expected pair surface density h(R) is given by 
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where n(z) is the cumulative quasar luminosity function down 
to a limiting magnitude (in this case i = 21), / CO mp(z) is the 
completeness in selecting binary candidates for follow-up 
spectroscopy as quantified in Paper I (see their fig. 7), and 
V c is the comoving volume between redshifts z m in and z m ax 
covered by the binary survey. The factor of 1/2 in eqn. (f2|i 
removes duplicate counts of pairs. 

Note we only consider the completeness in target selection 
/comp = /targ (i.e., the fraction of quasar binaries that would 
have been selected by the algorithm in Paper I) throughout 
this section. We discuss the effects of the completeness of the 
spectroscopic follow-up of survey candidates (i.e., the frac- 
tion of targets that have been observed), / spec , in ' 

Defining the usual quantity S = -2 In C we have 
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with all the model-independent additive terms removed. Here 
[^mi n ,^ max ] is the range of comoving scales over which 
we search for quasar pairs. To include all observed pairs 
with angular separation 9 < 60", we choose [R m in,Rnax] = 
[0.04,1] /i _1 Mpc for the low-z bin and [R m m,Rm^] = 
[0.1, 1.3] /r'Mpc for the high-z bin. We verified that our re- 
sults were not sensitive to the exact values of these limits. If 
we fit both ro,ML and 7ml we found that the best-fit model 
favors 7ml > 2.3 for both redshift bins, a value substantially 
steeper than the slope on large scales, 7^2 (e.g., Shen et al. 
2007). However, the spectroscopic completeness, f spec , prob- 
ably depends on angular separation, because we tended to ob- 
serve the closest candidates first; this may introduce an arti- 
ficially steep slope in the correlation function. Therefore we 
fix the slope 7ml = 2 (i.e., close to the measured slope of the 
large-scale correlation function, Shen et al. 2007) and mini- 
mize the merit function S with respect to ml only. A power- 
law slope 7 ~ 2 is also found for the clustering of SDSS LRGs 
to z — 0.4 (e.g., Masjedi et al. 2006) and photometric SDSS 
quasars (e.g., Myers et al. 2006, 2007a) over a wide range of 
scales down to r = 0.01 /i _1 Mpc . 

Alternatively, we may estimate the projected correlation 
function, i.e., the w p statistic, for these pairs. Following the 
definition of w p (e.g., Davis & Peebles 1983), we have 



w p (R) = 2 dH£ 2D {R,H) = 2 
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where tt(R\-R 2 ) is the projected comoving area of the cylin- 
drical annulus over which we search for pairs, R is the geo- 
metric mean pair separation in the bin, n(z) is the cumulative 
quasar number density, Af pa i r = YjDD is the observed num- 
ber of quasar pairs in the bin, and RR is the expected num- 
ber of random-random pairs in the cylindrical shell with radii 
(R\,Rz) and height dH. Note that there are some approxima- 
tions and ambiguities involved in Eqn. (|4j, such as the posi- 
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tion of the bin center, hence it can only be treated as a crude 



estimate for w p . 
For both the ML approach and the w p 
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FIG. 1 . — Measurements of the small-scale clustering for the low-z bin (left) and high-z bin (right). Filled circles are the large-scale correlation function data 
from Shen et al. (2007, the all sample) and dashed lines are their power-law fits with fixed slope 7 = 2. Squares are our estimate of w p using Eqn. |4), estimated 
in a large radial bin (filled) and two smaller radial bins (open). Points are placed at the logarithmic mean of pair separations in the bin, ho rizon tal error bars show 
the bin size, and vertical error bars show Poisson errors. The black hatched regions show our ML power-law fits to the small-scale pairs ( j|2.1l /spec = 1), with the 
vertical extent enclosing the 1<t statistical uncertainty from the M L fitti ng. If we assume minimal spectroscopic completeness, / spe c = 0.38 (0.52) for the low-z 
(high-z) bin, the ML results are shown as red hatched regions (see j|2.2t . These estimates, however, should be considered as solid upper limits. 

z bin, and ro,ML = 18. 22^ A ^ 2 /r'Mpc for the high-z bin, where 
errors are la statistical only; these results are shown as black 
hatched regions whose horizontal and vertical extent encloses 
the fitting range and statistical errors. For the binned w p statis- 
tic, we take all the pairs and use Eqn. (0 to estimate w p for 
the two redshift bins with Poisson errors. We then plot the w p 
estimates at the (geometric) mean values of separations (R) as 
filled squares in Fig. [T] To indicate the uncertainties in the bin 
center, we draw horizontal error bars which enclose the fitting 
ranges in the ML approach. In both redshift bins we further 
divide the pairs into two radial bins (with more or less equal 
number of pairs each), with the dividing scale R = 0.34 and 
0.56 /r'Mpc for the low-z and high-z cases (the dividing scale 
is set by the geometric mean of the maximum separation of 
observed pairs in the inner bin and the minimum separation 
of observed pairs in the outer bin). The w p estimates for the 
divided R bins are shown in open squares in Fig.Q] The results 
of the w p statistic are consistent with the ML results within the 
errors. However, due to the ambiguity of placing bin centers 
when there are only a few pairs, the w p data points cannot be 
used in the power-law fit. Our ML approach is not subject 
to such ambiguities, and therefore provides reliable clustering 
measurements. We tabulate the ML results in Table[T] 

2.2. Systematic Uncertainties 

Here we give some quantitative estimation of the systematic 
uncertainties in our ML results. The two major systematics 
come from the adopted luminosity function and the sample 
completeness. Our model luminosity function is quite uncer- 
tain down to i = 21, especially at z > 3.5 where no direct op- 
tical LF data are available. As we described above, the uncer- 
tainty in the LF is ~ 20%. In addition, the relative uncertainty 
in our pair target selection completeness is < 10% (Paper I). 
These taken together, introduce a systematic uncertainty in 
the best-fit r ,ML of o> = ±1.5/z _1 Mpc and ±3.1 /r'Mpc for 
the low-z and high-z bins, respectively; these values are com- 



statistic we need to 
This requires knowl- 
edge of the faint end of the luminosity function (z < 21) of 
quasars at redshift 2.9 < z < 4.5. We have searched the litera- 
ture for usable LF within these redshift and luminosity ranges 
(e.g., Wolf et al. 2003; Jiang et al. 2006; Richards et al. 2006; 
Hopkins et al. 2007). The Jiang et al. (2006) LF data probe 
sufficiently faint but do not extend to z > 3.6; the Richards 
et al. (2006) data have the desired redshift coverage but do 
not probe deep enough. By comparing the Richards et al. LF 
with the COMBO- 17 LF (Wolf et al. 2003), we found that 
the COMBO- 17 PDE fit gives better estimates of the LF at 
z > 3.5, e.g., it agrees well with the Richards et al. LF at 
the high luminosity end, and produces the expected flatten- 
ing at fainter luminosities. Motivated by these comparisons, 
we adopt a combination of the Jiang et al. fit (at z < 3.5) and 
the COMBO-17 PDE fit (at z > 3.5) for the model LF, scaled 
to our standard cosmology. We estimate an uncertainty in the 
cumulative number density (i < 21) of ~ 20%, based on the 
statistical uncertainties in these LF fits and comparison be- 
tween these optical LFs and the bolometric LF compiled by 
Hopkins et al. (2007), where the faint end LF at these redshifts 
is further constrained by X-ray data. This estimate of uncer- 
tainty in the model LF is, however, conservative at z > 3.5, 
since there are no direct optical LF measurements down to 
; = 21 within this redshift range. We will discuss the contri- 
bution of the uncertainty in the LF to the systematic errors in 
our small-scale clustering measurements, in 32.21 

Our clustering measurements are summarized in Fig. Q] 
where we plot for comparison the large-scale (R>2 /r'Mpc) 
correlation function data from Shen et al. (2007, the all sam- 
ple), for the low-z (left) and high-z (right) bins respectively. 
The ML approach yields to.ml = 8.3 1^} n 2 /i _I Mpc for the low- 
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TABLE 1 

Estimates of r for fixed power-law (7 = 

FUNCTIONS 



2) CORRELATION 





''O.ML (/spec = 1) 

/r'Mpc 


D,ML (lowest /spec) 
/r'Mpc 


ro (large-scale) 
/r'Mpc 


low-z 




i3.8i+| : |l 


14.79±2.12 


high-z 


18.22^ 


25 43 +4J6 


20.68 ±2.52 



FIG. 2. — HOD model predictions for a flux limit of ;' < 21 compared with the clustering data (notations are the same as Fig.[T}. The solid lines are the HOD 
predictions, where the dotted lines are the two-halo term contribution to w p . Three HOD models with satellite halo duty cycle f s = 0.1 (blue), 0.5 (cyan) and 1.0 
(magenta) are presented. For clarity, we have removed the upper limits on the small-scale clustering shown in Fig.[T](the red hatched regions). 

of the excess clustering occurs for the z < 3 sample (Hen- 
nawi et al. 2006), and 2) the quasar sample in Shen et al. 
(2007) has i < 20.2, while our binary sample has i < 21, thus 
luminosity-dependent clustering at such high redshift and lu- 
minosity ranges might play a role (e.g., Shen 2009). In the 
next section we show how these small-scale clustering mea- 
surements can be used to constrain halo occupation models. 

3. DISCUSSION 

The small-scale clustering measurements presented above 
can be used to constrain the statistical occupation of quasars 
within dark matter halos at z > 3. Given that we have a 
poor understanding of the physics of quasar formation, we 
use a simple phenomenological model relating quasars to ha- 
los to model the observed clustering results. The details of the 
model will be presented elsewhere (Shankar et al., in prepara- 
tion); below we briefly describe the model assumptions. 

We assume there is a monotonic relationship between 
quasar luminosity and the mass of the host dark matter halo 
(including subhalos), with a log-normal scatter S (in dex). 
Therefore for a flux-limited quasar sample, the minimal halo 
mass M m ; n and the average duty cycle /, defined as the frac- 
tion of halos that host a quasar above the luminosity threshold 
at a given time, can be jointly constrained from abundance 
matching and the clustering strength (e.g., Martini & Wein- 
berg 2001; Haiman & Hui 2001; Shen et al. 2007; White et 
al. 2008): 

r OO 

«qso,;<2i(z)= / /(M,z)$ ha i (M,z) 



NOTE. — The second column li sts ou r ML results assuming spec- 
troscopic completeness / spe c = 1 ( j)2. It . The t hird column lists ML 
upper limits assuming the lowest / spcc ( j|2.2> . The fourth column 
lists the large-scale correlation lengths from Shen et al. (2007, the 
all sample). Uncertainties are 1<t statistical only. 



parable to the statistical uncertainties reported above. 

In addition, our spectroscopy is incomplete even at 8 < 60" 
- only ^38% and ~ 52% of the high-priority low-z and high- 
z binary targets have been observed (see Table 3 of Paper I). 
Therefore we are undoubtedly missing some quasar pairs and 
our ML results are lower limits. Because targets further away 
from the stellar locus were assigned higher priority (Paper I), 
it is difficult to assess the effective spectroscopic complete- 
ness (because the most promising candidates were observed 
first); we expect that the effective spectroscopic completeness 
is larger than 50%. In the extreme case / spec = 0.3 8 (lo w-z) 
and 0.52 (high-z), we repeat our ML analysis in $2.11 with 
/comp = /targ x / spec and find r ,ML = 13.81^:^/! _1 Mpc and 
fu.ML = 25.43^2g/i _I Mpc for the low-z and high-z case re- 
spectively, where errors are la statistical. These estimates are 
shown as red hatched regions in Fig. Q] and should be consid- 
ered as solid upper li mits. 

The ML results in $2. II have comparable or lower cluster- 
ing amplitude at 0.1 < R < 1 /r'Mpc than the extrapolations 
from the fits for the large-scale correlation functions (Shen et 
al. 2007, 2009). This does not directly contradict the results 
in Hennawi et al. (2006) for z < 3 quasars since: 1) our sam- 
ple barely probes scales below R ~ 0.1 /r'Mpc where most 
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where «qso,;<2i(z) is the cumulative quasar number density 
with flux limit i < 21, M is the halo mass, $haio(M,z) is the 
halo mass function per logM interval, and < /(M,z) < 1 is 
the average halo duty cycle, which may be a function of both 
redshift and halo mass. 
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In general the halo mass function <&ha\o(M,z) includes con- 
tributions from both halos ($ c ) and their subhalos ($ s ), where 
we use the Sheth & Tormen (1999) halo mass function for the 
former and the unevolved subhalo mass function from Gio- 
coli et al. (2008) for the latter. It is important to use the 
unevolved mass (i.e., mass defined at accretion before tidal 
stripping takes place) for subhalos, since subhalos will lose a 
substantial fraction of mass during the orbital evolution within 
the parent halo. We denote the average duty cycles for central 
and satellite halos as f c and f s respectively. Note that we as- 
sume halos and subhalos of the same mass host quasars of the 
same luminosity - of course, subhalos within a given halo will 
be less massive and thus host quasars fainter on average than 
the central quasar. The satellite duty cycle /, is the fraction 
of black holes in subhalos that are active at a given time. The 
fraction of luminous quasars that are satellites is always small, 
regardless of f s , because the number of massive satellite halos 
is itself small. 

An important consequence of the rareness of binary quasars 
is that the abundance matching, i.e., Eqn. (O, can be done 
using central halos only, and we have / w f c , <l>h a io ~ 3>c m 
Eqn. (0; the satellite duty cycle f s will only affect the small- 
scale clustering strength. In order to simultaneously match the 
large-scale clustering of z > 3 quasars (Shen et al. 2007, 2009) 
and their abundance, Shankar et al. (2008, 2009) found large 
values of duty cycle f c ~ 0.5 - 1 are needed, as well as small 
scatter for the quasar-halo correspondence, if the Sheth et al. 
(2001) bias formula is used (cf. Shen et al. 2007 for the usage 
of alternative bias formulae). For simplicity we fix f c = 1 and 
£ = 0.03 dex in what follows, which produces adequate fits 
for the large-scale clustering and abundance matching 13 . 

Using Eqn. (O, we determine the minimal halo mass to be 
M m i n ~ 1O 13 /! _1 M for both redshift bins. We then use stan- 
dard halo occupation distribution (HOD) models (e.g., Tinker 
et al. 2005) to compute the one-halo term correlation function 
with different values of satellite duty cycle < f s < 1 . 

Fig. [2] shows several examples of our HOD model at z = 3.1 
(left panel) and z = 4 (right panel) with /, = 0.1 (blue), 0.5 
(cyan) and 1 .0 (magenta) for a flux limit of i = 2l. Solid lines 
are the total correlation while the dotted line denotes the two- 
halo term contribution. As expected, the value of f s has no 
effect on the large-scale clustering; it only changes the small- 
scale clustering amplitude. These are not actual fits to the 
data because the quality of our measurements does not allow 
a reliable HOD fit. Nevertheless, it seems that some active 
satellite halos are required, but only < 50% of satellite ha- 
los can be active at a given time in order not to overshoot the 

13 Although the model still underpredicts the large-scale clustering a bit 
for the high-z bin even with f c = 1, as noted in earlier papers (White et al. 



small-scale clustering. This constraint is less stringent if we 
consider instead the upper limits on the small-scale clustering 
discussed in 32.21 One potential concern regarding our model 
is that the adopted subhalo mass function has not yet been 
tested against simulations for the extreme high-mass end and 
redshift ranges considered here; nevertheless our model ap- 
proach demonstrates how the small-scale clustering measure- 
ments can be used to constrain quasar occupations within ha- 
los. We defer a more detailed investigation on the uncertain- 
ties and caveats of our halo models to a future paper (Shankar 
et al., in preparation). 

4. CONCLUSIONS 

We have measured the small-scale (0.1/r'Mpc < R < 
1 /r'Mpc) clustering of quasars at high redshift (z > 3), based 
on a sample of 15 close binaries from Paper I. Strong clus- 
tering signals are detected, comparable to or lower than 
the extrapolations from the large-scale clustering based on 
SDSS quasar samples. The small-scale clustering increases 
in strength from z ~ 3 to z ~ 4, consistent with that of the 
large-scale clustering (Shen et al. 2007, 2009). 

Using a simple prescription relating quasars to dark matter 
halos, we constrain the average duty cycles of satellite ha- 
los at z > 3 from the small-scale clustering measurements. 
We found tentative evidence that only ~ 10% -50% of satel- 
lite halos with mass > 10 /rtif© can host an active quasar 
(with i < 21). With the completion of our ongoing binary 
quasar survey, we will have better estimates of the spectro- 
scopic completeness and therefore will confirm our results. 

Future surveys of fainter binary quasars at z > 3 will 
increase the sample size and hence the signal-to-noise ra- 
tio of the small-scale clustering measurements. These 
measurements, together with better understandings of the 
halo/subhalo abundance and clustering at z > 3 from simula- 
tions, will provide important clues to the formation of quasars 
at high redshift. 
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